% SIMPLE 2-D PROBLEM
numel = 1;
Y = 10;
it = 1;
if it == 0
    q = 1;
else
    q = 1;
end

dat = [ 1.500   0.030   0.000  1.0 0.000  0.000  
        0.420   0.080   0.125  0.0 0.020  0.000
        2.000   0.040   0.000  0.0 0.000  0.000
        0.300   0.010   0.000  0.0 0.040  0.000 ];     
abdat = [0.010
         0.080
         0.000
         0.010];    

numg = 2;  numm = 2;

if numel == 1
% still simple, but nine coarse meshes, uniform
xcm  = [ 0   2   4   6 ];
xfm  = [   1   1   1   ]*10;
ycm  = [   0     1  ];
yfm  = [     3      ];
mt   = [   1   2   1   ]';  
else
%still simple, but nine coarse meshes, uniform
xcm  = [   0     2     4      6  ];
xfm  = [     1      1      1     ]*Y;
ycm  = [   0     2     4      6      8     10    12];
yfm  = [     1      1      1     1      1     1    ]*Y;
mt   = [  1 1 1 1 1 1
          1 2 1 1 2 1
          1 1 1 1 1 1 ];         
end

BCL = 0;
BCR = 0;
BCB = 1;
BCT = 1;        
       

IBSL = [ 0 
         0 ]';
IBSR = [ 0     
         0 ]';  
IBSB = [ 0                   
         0 ]'; 
IBST = [ 0
         0 ]';     

src=0;